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Abstract 



The electron neutrino mass has been measured in several tritium beta de- 
cay experiments. These experiments are sensitive to a small neutrino mass 



> 

o 

■ because the energy release of the decay is small. But the very smallness of 



the energy release implies that the Coulomb interactions of the slowly moving 
emitted beta electron are relatively large. Using field theoretic techniques, we 
derive a systematic and controlled expansion which accounts for the Coulomb 
^ . effects, including the mutual interaction of the beta ray electron and the elec- 



tron in the final ^He"*" ion. In our formulation, an effective potential which 
describes the long range Coulomb force experienced by the beta ray is intro- 
duced to ensure that our expansion is free of infrared divergences. Both the 
exclusive differential decay rate to a specific final ^He"'' state and the inclu- 
sive differential decay rate are calculated to order rf^ where rj is the usual 
Coulomb parameter. We analyze the order rj^ correction to the beta ray spec- 
trum and estimate how it may affect the neutrino mass squared parameter 
and the endpoint energy when this corrected spectrum is used to compare 
with the experiments. We find that the effect is small. 



*Present address: Department of Physics, Purdue University, West Lafayette, Indiana ^7907 

1 



I. INTRODUCTION AND SUMMARY 



The mass of the electron's anti-neutrino m,^ has been measured in several tritium beta 
decay experiments . These experiments are quite sensitive to a possible neutrino 

mass since the end-point energy Q in this beta decay is very small {Q ^ 18.6 KeV), and 
an examination of the beta ray energy spectra near the end point region can thus reveal 
a small neutrino mass. On the other hand, the very smallness of the decay energy implies 
that the Coulomb interactions between the outgoing beta electron and the atomic electron 
in the final helium ion may play a significant role in the interpretation of the experiments. 
In units where c = 1 and h = 1, the strength of this Coulomb interactions is governed by the 
parameter 



where a ~ 1/137 is the fine structure constant and v is the emitted beta electron velocity, or, 
equivalently, where Oq is the Bohr radius and p is the momentum of the outgoing beta ray. 
Since the end-point energy Q is small, the beta electron velocity is small and the Coulomb 
parameter rj is relatively large. Near the tritium decay end point, rj ^ 0.03. 

To study a small neutrino mass, modern experiments use data from the beta ray spec- 
trum end point to approximately 800 — 1000 eV below the end point The actual 
experiments work with molecular tritium. The theoretical formula for the spectrum used in 
the analysis of the experimental data is the sudden approximation result for the molecular 
final states which neglects the Coulomb interaction between the beta ray electron and the 
electrons in the daughter '^He — T+ molecule 0,H,§,|^,PJ^,|] . In this paper, we shall investigate 
the Coulomb interaction effect, hereafter referred to as the atomic effect, and estimate its 
size. Instead of molecular tritium, we calculate the effect for atomic tritium, where the cal- 
culation is simpler and brings out the essential points. The calculation for molecular tritium 
would be parallel to that of atomic tritium, but much more complicated. Currently, tritium 
beta decay experiments appear to find a squared electron neutrino mass which is negative. 
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ml ^ — 60eV^ . It should be noted that there are two relevant dimensional pa- 

rameters which enter here: The atomic energy scale Ry ~ 13.6 eV and the energy range AE 
where the experimentalists fit their data. Typically, AE ~ 10^ eV. Thus, apart from po- 
tentially large numerical factors, the effect of the atomic corrections on the determination of 
the the neutrino mass squared may be of order rfAE"^ ^ 900eV^ as well as rf^"^ — 0.2eV^. 
We shall provide a detailed analysis of how much the atomic effect influences the neutrino 
mass squared parameter. 

The sudden approximation spectrum used for fitting the experimental data takes the 
form 



AFi2,E)p Y: WfUQf-E)^iQf-E)^-ml, (1.2) 

/</m ax 

where E and p are the energy and momentum of the outgoing beta electron and ml is the 



neutrino mass squared. The spectrum ( |1.2| ) needs some explanation. The first factor A is 
an overall constant, and F{2, E) is the usual Fermi function, which is the square of the beta 
ray wave function evaluated at the ^He"^ nucleus, for the Z = 2 charged helium nucleus. The 
atomic transition probability Wfi for the initial tritium state \i) decaying to the particular 
final ionic state (/| is the squared matrix element 

Wf, = \Tf,\' ^ (1.3) 

The energy release 

Qf = Q + E,-Ef (1.4) 

is that for the decay process which starts with the initial tritium atomic energy Ei and ends 
with the final ^He^ ion energy Ef. The upper limit /max on the sum over states is set by the 
conservation of energy. It mostlyB corresponds to a state of '^He^ that has become unbound 



^In a small region at tfie very end of tfie beta ray spectrum, this upper limit corresponds to a bound 
state. We ignore these cases since statistically this small region makes only a small contribution 
which does not significantly alter our results. 



into ^He^"*" + e with the e kinetic energy being 

E^^^ = K^/2m = Q + Ei-E. (1.5) 

Parameters including the overall strength factor A, the end point energy Q, and the neutrino 
mass squared ml are determined by fitting the spectrum ( |1.2| ) to the measured spectrum. 

To simplify our work, we shall neglect relativistic corrections and treat the beta ray non- 
relativistically. This is valid since relativistic effects, of order f ^ ^ 0.06, provide only small 
corrections that can easily be accounted for. They correct both the sudden approximation 
result and the atomic effect. The correction to the sudden approximation result may be 
taken care of, to good accuracy, by using an approximate relativistic Fermi factor F{Z, E) in 
Eq. ( |1.2D The correction to the atomic effect is of the order r/^f^, which is negligible 



In calculating the beta ray decay rate, the usual perturbative expansion [|Tl|,|l3 uses the 
Coulomb potential produced by the final helium nucleus as the zeroth-order approximation 
for the potential experienced by the emitted beta electron. This zeroth-order approximation 
does not describe the long distance behavior of the potential in which the real beta ray 
moves since the screening effect of the outside electron has not been taken care of. Therefore, 
infrared divergences appear and make the usual perturbative calculation inconvenient. Using 
field-theoretic methods involving the reduction technique, we shall instead make use of a 
conveniently chosen effective potential as the zeroth-order approximation. Since this effective 
potential properly accounts for the long-range character of the screened Coulomb potential, 
we can perform a systematic expansion in powers of the small parameter rj, with the expansion 
free of infrared divergences. We shall compute the first term in this expansion, which is of 
order 77^, for the decays into the individual final ionic states. The inclusive sum of these 
individual decay rates agrees with previous results [pT| , |l2| under two approximations! made 



^These approximations are made only for the order r/^ correction to the exclusive decay rates 
caused by the atomic effect. 
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there: 1) a "uniform phase space factor approximation" and 2) the closure approximation. 
The uniform phase space factor approximation takes the neutrino phase space factor, which 
depends on the final helium ionic energy, to be the same for different final helium ionic states. 
The energy of the final ionic state appearing in the phase space factor is replaced by the 
helium ionic ground state energy [|TT|. Within this uniform phase space factor approximation. 



references [|ri|,0 then make the closure approximation, which extends the upper limit /max 
of the sum over final helium ionic states to infinite energy. The closure approximation over 
counts contributions from the final helium ionic states whose energies lie above the limit set 
by energy conservation. In most regions of the beta ray spectrum, this approximation is valid 
because the limit determined by energy conservation is much higher than the atomic energy 
scale. However, the closure approximation cannot be used for the region near the end point of 
the beta ray spectrum since here little energy remains to excite the final helium ionic states, 
and threshold effects must be taken into account. Since the shape of the beta ray energy 
spectra near the end point region is crucial for the neutrino mass measurement, we shall not 
rely on either the closure approximation or the uniform phase space factor approximation. 

Our result for the case of a vanishing neutrino mass, m,^ = 0, may be expressed in the 
form 

n\- 777 r -1 

-^ = ^,F{2,E)p\Tf,\\l-2v')[P{Q-E)+v'R{Q-E)] , (1.6) 
where m is the electron mass, the Fermi factor 

^(2,^) = Y^Pi^ (1-7) 

is the value of the squared non-relativistic beta electron wave function at the ^He"*" nucleus, 
and |T^P is squared nuclear decay amplitude in the absence of atomic corrections. Here and 
in the following we shall make use of the Rydberg energy 

IV = -a^ m = — ~ 13.6 eV . (1.8) 
2 2ao ^ ' 

In the closure approximation which we have just discussed. 



T.\TM' = l, (1.9) 

all / 

and a simple calculation, reviewed in Section IV, gives 

Y: -Er = {Q~E + 2Byr + 4 {Byf . (1.10) 

all / 

This provides an approximate evaluation of the sum in Eq. (|1.2|) when = 0. The finite 
upper limit on the sum alters this evaluation and defines the first function that appears in 
the curly braces in Eq. (|1.6|), 



P{Q-E)^ Y: iTf^Qf-Ef. (1.11) 

/</m ax 

For later convenience, we take the independent variable to be Q — E, which is related to 
-E-max = K^/^m by Eq. ( p..5|) . In Section IV we compute the correction S{Q—E) to the 
closure approximation such that 

P{Q-E) = {Q-E+2 Byf+4: By^+S{Q-E) . (1.12) 

The term 7fR{Q—E) in Eq. (|1.6|) represents the correction due to the atomic effect. The 
major purpose of our work is to compute R{Q—E). 

In order to describe the effect of our result, we shall slightly simplify the formula (|1.2| ) 
used in the experimental data analysis. Although the experiments measure the spectra near 
the end point E = Q, they nonetheless mainly measure electron energies with {E — QY ^ 
since the decay rate is so small in the region very near the end point. Hence one may expand 
the square root in Eq. ( |1.2| ) to get a simpler form: 



dV 



AF{2,E)p Y: \Tfr? {Qf-Ef--m 



dE 



1 



,2 



2 

The first sum which appears here is just the sum (|1 . 1 1|) , and defining 



:i.i3) 



Pi(g-E)= Y \Tf^\\ (1.14) 

/</m ax 

which is the total probability for the initial tritium to make the transitions to any final state 
below energy -Emax, we have 
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^ = AF{2,E)p 



P{Q-E) - ]^nilPr{Q-E) 



:i.i5) 



To see how the atomic correction affects the neutrino mass squared, we note that Eq. (|1.6| ) 
reduces to the sudden approximation result (|1.15|) for ml = if R{Q—E) is ignored. How- 
ever, if we include the atomic correction R{Q—E), and require that formula (|1.2| ) fits data 
described by spectrum (|1.6|) with free parameters A, Q, and m^, we may get a nonzero mf, 
for the best fit. The correction 7fR{Q—E) to the spectrum effectively changes A, Q, and 
ml by AA, AQ, and Am^ in such a way that the change 



A A 1 

P{Q-E) + AQ P'{Q-E)--Aml P,{Q - E) 



:i.i6) 



to Eq. (|1.15| ) mimics the correction to the spectrum represented by rf'R{Q—E). Here P'{Q — 



E) is the first derivative of the function P, and we have ignored terms higher order in the 
small parameters AA, AQ, Ami, m^. To find the changes Ami ^^^1 AQ due to our 
atomic effects of order r]"^, we fit the atomic correction represented by the functioni R{Q—E) 
to a linear combination of P{Q — E), P'{Q — E), and Pi{Q — E) as in the form ( |1.16| ). The 
neutrino mass measurement is sensitive to energies from the beta ray spectrum end point 
down to approximately 59 — 74 (800 — 1000 eV) below the beta ray spectrum end point 
P,0,^,||J^. Fitting in this energy range, we find in Section VI that 

R{Q-E) ^ OM P{Q - E) - A.6By P'iQ - E) + 6.3Ihf^ Pi{Q - E) , (1.17) 

which, in view of Eqs. ( |1.6[ ) and ( [l.l(j| ), gives 



Am^ ^ -12.6r/2%2 ~ -1.7eV% (1.18) 

and 



^Although T]'^ depends on the beta ray energy, we shall treat it as a constant in the region near 
the endpoint used to determine the neutrino mass. This is valid because rj"^ varies only by a small 
amount ~ rj^AE/Q as the beta ray energy varies from the end point to AE below the end point. 
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Table 1: The dependence of the neutrino mass squared fitting on the range of the energy 
where the fit is made. 





range 

o 


of -Emax 


Ami 


- 


- 49 IV 


= - 


- 670 eV 


- 0.80 eV^ 


- 


- 55 IV 


= - 


- 740 eV 


- 1.2 eV^ 


- 


- 59 % 


= - 


- 810 eV 


- 1.5 eV^ 


- 


- 64 % 


= - 


- 870 eV 


- 1.7 eV^ 


- 


- 69 % 


= - 


- 940 eV 


- 2.1 eV^ 


- 


-74% = 


= — 


- 1010 eV 


- 2.5 eV^ 


- 


- 79 % = 


= — 


- 1080 eV 


- 2.8 eV^ 



AQ ^ -4.6 ri^'Ryc^ -0.047 eV . 



:i.i9) 



The fitting formula ( |1.17D depends on the range of energy where we do the fit. Equation (|1.17| ) 
is done with the energy range -Emax being — 64 Ry (0—870 eV). Increasing the range of energy 
for Einax from 0-49% (0-670 eV) to - 79% (0-1080 eV) changes the parameter Ami 
linearly. For each 5 % increase of the energy range, Ami decreases by 0.4 eV^. We show this 
dependence in Table 1. The atomic effect causes only a few eV^ correction to the neutrino 
mass squared and thus does not affect current experimental bounds on the neutrino mass 
squared. 

The next five sections describe our methods and calculations. In Section II, we develop, 
using quantum field theory techniques, a general formula which is convenient for calculating 
the exclusive decay rate to a specific final '^He''' ion state. A comparison potential is intro- 
duced to facilitate a systematic and controlled expansion in the small parameter rj. With 
an appropriate choice of the comparison potential, given in Section III, the exclusive decay 
rates are calculated to rj'^ order in Sections IV and V. Finally, in Section VI, we examine the 
atomic effects in the determination of the neutrino mass squared. Fitting formulas for the 
beta ray spectrum are given for the region near the end point. In Appendix A, we calculate 
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the value of the wave function of the beta ray at the origin for the sudden approximation. 
Our evaluation of this amplitude confirms an old result [|1^] and shows, incidentally, that it is 
accurate through terms of order rf rather than being valid only through order rf. The details 
of evaluating the second order atomic correction are presented in Appendix B. In Appendix 
C, we provide the details of the computations for summing over the final states. Finally, 
in Appendix D, we investigate the nature of the exchange corrections. These correspond 
to the process in which the electron produced by the weak interaction is bound in the final 
^He^ ion with the initial tritium atomic electron being ejected. We show that the leading 
exchange corrections to the amplitude are of order yf, in contradiction with the previous 



order rj^ estimate that has appeared in the literature [T^. The exchange corrections to the 
decay rate, however, are of order 77^. 

We should emphasize that, although the corrections which we have found for the neu- 
trino mass determination from tritium beta decay are not significant, the methods which we 
have developed to treat the problem may be useful in the examination of other beta decay 
processes. 



II. GENERAL FORMULATION 

Since the emitted beta ray has a maximum energy of 18.6 KeV, in natural units where 
c = 1, the square of the electron velocity is less than 0.1, it is valid to treat the electron non- 
relativistically. The interactions due to the spins of electrons and nucleons are relativistic 
effects. Therefore, the spectrum is not significantly affected by the spins of these non- 
relativistic particles except for an overall constant which accounts for the spin degrees of 
freedom. Hence, the effective Hamiltonian for the beta decay for our purpose may be written 
in the form 

Hf^ = (const) I (dh) ^t(r)^„(r)^t(r)^.(r) + h.c. , (2.1) 

where the operators ipp, i^n, i'e, and ipf are spinless field operators, which destroy the proton, 
neutron, electron, and neutrino, respectively. 
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We will calculate the decay rate of the tritium atom expressed as an energy spectrum 
integral of the beta ray. To accomplish this, we assume the following process: Initially, at 
time t = Ti, there is a tritium atom in its atomic ground state denoted by at a later time 
t — —T/2 the weak interaction is turned on and then turned off at time t — T/2; finally the 
outgoing electron wave packet travels for certain time and hits the detector far away from 
the tritium atom at time t = Tf. Here T is the time for the interaction to act, and the 
relations Ti < —T/2, T/2 < Tf are assumed. We shall then calculate the probability for 
detecting the outgoing electron, or the decay rate by dividing this probabihty by T. 

To first order in the weak coupling, the amplitude for the tritium atom decaying to a 
specific final state is 

rT/2 

AT^-t{u,eJ;Tf\ j_^^^dt H^^{t) \i;T,) , (2.2) 

where we have denoted the initial tritium atom state at time t = Tj by \i;Ti) and the final 
state which contains an anti-neutrino, an outgoing electron, and the helium ion ^He"^ at time 
t — Tf hj (z/, e, f;Tf \ . The time dependence of the Hamiltonian H^(t) is governed by the 
Hamiltonian H which contains the particles kinetic energies and the Coulomb interaction 
between the charged particles. 

In the infinitely heavy nucleus limit, which is appropriate in this non-relativistic situation, 
the operator ■0^(r, t)'0„(r, t), when evaluated between states (/| and is proportional to 
5(r)(5+(t), where Q+{t) is the charge raising operator which converts a neutron into a proton, 
and the origin is chosen to be located at the nucleus. Under this infinitely heavy nucleus 
approximation, the amplitude for the tritium beta decay reduces to 

/•T/2 

At = (const)(-^) j_^^^dte'''^'{e, f-Tf\4{0,t)Q+{t)\t-T,) , (2.3) 

where the spatial integral has been performed by exploiting the 5-function produced from 
the infinitely heavy nucleus limit, and the action of the neutrino field operator in the weak 
decay effective Hamiltonian has been used, 

(P, e, /; TflMO, t) = e^^''*(e, f;Tf\. (2.4) 
10 



Since only the electron non-relativistic field operator is involved in the following calculations, 
we shall change notation and replace ipg hy ip. 



A. Reduction Technique with a Comparison Potential 

We now apply the reduction method to the final outgoing electron. To achieve this, we 
first construct an asymptotic wave packet for the final outgoing electron, which propagates 
in the Coulomb potential of the produced '^He^. Since the Coulomb potential is a long range 
interaction, the wave packet of the final outgoing electron must be constructed as moving 
in the Coulomb potential or a comparison potential v{r) having the same large distance 
behavioiS, 

v(r) ^ —- — as r — > oo . (2.5) 
47rr 

Since we only require that the comparison potential v{r) has the correct long distance be- 
havior, the result ( p.26| ) which we shall obtain should not depend on the specific choice of 



v{r). This we shall prove later. In view of these considerations, the final electron-ion state 
can be expressed by the asymptotic state 

(e, /; Tf\ = (/; T;| | J (rf^r)$;^(r, T^)^(r, T;)} , (2.6) 

where the wave function (r, t) satisfies the Schrodinger equation 

"4^-^'''^ = {-^ + ^W} ■ (2-7) 

This electron wave packet can be viewed as the superposition 



^For the case where the final ^He"*" is in an unbound state, besides the anti-neutrino, the final 
state contains '^He''"^ and the two unbound electrons. The electron moving faster is identified as 
the beta ray electron. Therefore, the other electron which moves relatively slowly still plays the 
role of screening the ^He"'"'' nucleus. 
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$;(r,t) = l{d'p)gip)cp;iv)e'^\ (2.8) 



where g{p) is a probability amplitude peaked at p = Pe with a width Ap, and 0p(r) obeys 



2 



|-^ + ^(O|0;(r) = i5 0;(r), (2.9) 

with E = p^/2m and 0p(r) satisfying the boundary condition that for large r its asymptotic 
formi contains a plane wave with the momentum p. 

Since the time dependences of the final ^He"*" state and the initial tritium atom state are 
irrelevant phase factors, from now on we shall replace {f',Tf \ by (/| and \i;Ti) by Using 
the construction ( |2.6| ), the usual reduction method [|1^] gives 

,. ... d 

-T/2 



{d'^)j^ dx'j^^dte''^^' — {f\T{%^{x)ij{x)^\Q,t)Q+{t))\i), (2.10) 



where we have used 



II (rf3r)$;^(r,T,)V^(r,r,)} = , (2.11) 



since, at the very early time t = Ti, the wave function represents an incoming electron 
wave-packet which has no overlap with the initial localized electron state in the tritium atom. 

Using the superposition ( p.8|) and shifting the time integral variable to x° — t in 
Eq. (|2.10|) , the Heisenberg equation of motion then displays the t-dependence in the form 
exp{—iHt} ■ ■ ■ exp{iHt}, with the left exponential factor acting directly on the final state 
and the right factor acting directly on the initial state. Hence 

At = (const) (-0 Jid^x) j [d^v] 9{v) J^^\t e''~^f+^+^^-^^-'^^' 
d 



^;(x)e^^^' 



(/|T(v^(x)^t(o,o)g+(0))K)] , (2.12) 



dx^ 

where Q is the total energy released which equals the tritium-helium nucleus mass difference 
minus the electron mass, Ef is the atomic energy eigenvalue of the produced ^He^ ion, and 
Ei is the energy of the initial '^H atom. The time integration now yields 



^The asymptotic form of this energy eigenstate also contains an incoming wave which is the time 
reversal of the scattered wave. 
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At = (const)(-i) j (d=^p) ^(p) 27: 5T{Ef+E+E,-Ei-Q)r{p) , (2.13) 
where the function 5t{x) is defined by 

TTX 

which is approximately a "5-function" but with a width 1/T, and the atomic matrix element 

is defined as 

r(p) = J (d^x) A {0;(x)e^^^"(/|T (^(a;)^t(o)g^(o)) |^)} . (2.15) 

To calculate the atomic matrix element T[p), we write the time-ordered product in terms 
of step functions and take the partial derivative with respect to to get 

^(P) = <P;mf\Q+\^) + J (d'x) (/|T A (0;(x)e^^-V(x)) V'^(0)g+(0)|z) , (2.16) 

where we have used the equal-time anticommutator of the electron creation and annihilation 
operators, 

{V'(r),V'^(0)} = 5(r). (2.17) 

The time evolution of the electron destruction operator ip is governed by the Heisenberg 
equation of motion 

t^Mr, t) = l-^ + V{r, t)U(r, t) , (2.18) 



dt 2m 



where 



with 



p{r',t)=i^^{r',t)i^{r',t), (2.20) 
and the charge operator Z having the properties 
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ZQ+ = 2Q+, Q+Z = Q+, (2.21) 

since the charge of the tritium nucleus is unity and Q+ is the charge raising operator. 

Exploiting the equation of motion ( p.l8| ) of the operator ip^x) and the differential equa- 
tion (|2.9| ) satisfied by 0p(x) and integrating by parts yield^ 

T{p)=cf>;{0){f\Q^\^)-^^'dt^{dh){f\^c|);{r^^^^^ 

(2.22) 

Defining an operator Np(t) by 

iVp(t) ^ / (dh) 0;(r) [Vir, t) - v{r)]i^{v, t) , (2.23) 



which accounts for the difference between the real interaction experienced by the beta ray 
and the comparison potential, and utilizing 

Np{t) = e'"'Np{0)e-'^' , (2.24) 

we can do the time integral and obtain 

1 _ pi{E-Ei+H)Ti 

+ (/|^t(o)g^(0) iVp(O)K) . (2.25) 

Since we require that the potentials V{y) and v{r) have the same asymptotic behavior at 
large distance, the operator A^p(O) is localized near the origin. The facts that the final state 
(/I is localized and that the operator A^p(O) has a compact support, enables us to drop the 
highly oscillating terms in Eq. ( ^.25| ) .B Thus we obtain 



^The surface integral can be omitted since when the atomic matrix element T(p) is inserted into 
Eq. ( p.l3| ) to evaluate the amplitude At, the integral over the momentum p makes the integrand 
vanish on the surface by the Riemann-Lebesgue lemma if the surface is sufficiently large. 

gain, if we put the expression ( |2.25| ) for the atomic matrix element in Eq. ( |2.13| ), by the Riemann- 
Lebesgue lemma, the integral over momentum p makes the highly oscillating terms vanish. 
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r(p) = 0;(o)(/|Q+K) - (/|iVp(o)-^-^-^— -^t(o)Q^(o)|^) 

^f\^\0)Q^{0)j^^-^^^-—N^mz) , (2.26) 

where e is an infinitesimal positive number.! This is our general formula for calculating the 
atomic matrix element T(p) which exploits the comparison potential. For later reference, we 
shall denote the three terms in Eq. ( |2.26| ) by Ti, T2, and T3, so that this result is expressed 

as 

r(p) =Ti + T2 + T3. (2.27) 



B. Comparison Potential Invariance Theorem 



We now prove explicitly that the atomic matrix element T{p) does not depend on the 
short range behavior of the comparison potential v{r). To accomplish this, we consider a 



small variation Sv{r) in the comparison potential and show that T(p) evaluated by Eq. ( 2.26 ) 
is unchanged under this small variation. Varying v{r) in Eq. ( p.9|) yields 



p ■ 



Inserting the variation ( p.28| ) into the definition ( |2.23| ) of A^p gives 

^2 



2m 

2 

V 



^(r) 



J{dh)6<p;{r)lv{r)-^-EU{r 



(2.28) 



(2.29) 



This e specifies the correct way to avoid the zero in the denominator. Before dropping the 
highly oscillating term, while doing the p integral, the contour can go on the top of the zero of 
the denominator since the numerator is also zero and thus the integrand is finite. To apply the 
Riemann-Lebesgue lemma, we can deform the contour in the p plane slightly to avoid the zero in 
the denominator. As Tf ^ 00 and Tj —00, we can drop the highly oscillating terms with the 
contour properly deformed so as to make these large time limits well defined. 
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where we have performed two partial integrations and dropped the surface terms to make 
the differential operator V acting to the right. The dropped surface terms vanish because 
6v{r) must be a localized quantity so as to keep the correct long distance behavior in the 
comparison potential v{r). Noticing that 

< -^ + F(r)|^(r) = [^,i7], (2.30) 

we have 

5iVp(0) = \^Ji£r)5<p;{r)^Piv),H^-Elidh)S<p;{r)^Piv). (2.31) 

Utilizing this variation of the operator A^p(O), we indeed find a vanishing variation of the 
atomic amplitude, 

6T{p) = 6<p;{0){f\Q^\z) - (/I J {dh)6<p;{r)^{r)^\0)Q40)\t) 
-(/|V^t(o)g^(0) J{d'r)d<p;{T)^{r)\z) 
= 0, (2.32) 



upon using the anticommutation relation ( |2.17| ). 



C. Decay Rate 

We now take the plane wave limit to compute the decay rate. Taking the limit where 
the distribution amplitude g{p) becomes a delta function, completing the p integral in the 
expression (|2.13|) , and summing over all the possible final electron states expresses the decay 
rate as 

^ f (rf3p^) I- (rf3p) i(const) 27r5TiEf+E+E,-E,-Q) r(p)p 

^ " J (27r)3 J (27r)3 T " ^^"^^^ 



Noting that 



j.^ 2^ MEf+E+K-E,-Q)r ^ s{Ef+E+K-E,-Q) (2.34) 
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produces the decay rate 

r = /^/|^2^5(i5;/+i?+i?.-i?.-g)|r,nT(p)r. (2.35) 

Here we have compensated for our previous neglect of the real structure of the beta decay 
interaction by inserting a factor of the true squared amplitude |T^P in the absence of the 
atomic effects.! (In the absence of atomic effects, T(p) ^1.) Performing the momentum 
integrals yields the differential rate for the decay to a specific final ion state. Assuming that 
the neutrino mass vanishes, this gives 

dTfi m 



p{Qf-EY\Tp\'\r{jp)\\ (2.36) 



dE 27r3 

where m is the electron mass, and Q/ is the effective reaction energy for a specific final 
atomic state (/| defined by 

Qf = Q + E,-Ef. (2.37) 

When the decay takes place to an excited ^He^ ion with non-vanishing angular momentum, 
one must average over the spatial orientation of this state; this average is indicated by writing 



III. SECOND ORDER EVALUATION 



We turn now to study each term in our new expression ( p.27|) for the atomic amplitude 



and estimate their order in terms of the Coulomb parameter rj defined by 



^Summing over the final polarization of '^He nucleus, this squared amplitude is simply 

\Tfs\^ = GIcos^9,{\Gv\^ + 3\Ga\^), 

where Be is the Cabibbo angle and Gy — 1 and Ga — 1-25 are two coupling constants for the 
charged vector-current and axial-vector-current nuclear matrix elements which are nearly equal to 
the corresponding values for free neutron decay. 
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where Cq is the Bohr radius and a is the fine structure constant. The first term, 

Ti = 0;(o)(/|g+K) = 0;(o)T;,, (3.2) 

which is the result of the sudden approximation, dominates the amphtude. This term is of 
order 1. It can be calculated by explicitly choosing a comparison potential v{r) and solving 
an eigenvalue problem. 

The second term, T2, cannot be calculated analytically without doing any approximation. 
We shall find an appropriate approximation to calculate it for the case where the final 
outgoing electron has an energy much larger than the typical atomic energy scale so that 
the Coulomb parameter rj is small. We shall, in fact, compute the atomic amplitude through 
order 77^. It is convenient to rewrite the field theoretic expression of T2 in Eq. ( |2.26| ) in terms 
of ordinary quantum mechanics notation. Using this notation, the two electrons are labeled 
by the subscripts 1 and 2 and the action of the field operators produces states that are 
explicitly antisymmetrized. For the bra and ket, the first quantum number refers to electron 
1 while the second one refers to electron 2. With this notation in hand, T2 may be expressed 
as 

2e2 



T2 = -(p,/| 
where 



47r|ri — Attvi 



viri 



H - E - Ef -ie 



|ri=0,z) - K,r2=0)] , (3.3) 



(P,/I = / (d'v,) I (rf3r2)0;(ri)0}(r2)(ri,r2| , |ri=0,^) = J {dh2)\r,=0,r2)MT2) , (3.4) 



and |i,r2=0) is simply |ri=0, z) with 1 and 2 exchanged. Equation (|3.3|) naturally defines 
the direct term as the part associated with the ket |ri=0,i) and the change term Tg as 
the part associated with \i, r2=0) so that 

T2 = Td + Te . (3.5) 
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The Hamiltonian in Eq. ( p.3| ) may be obtained by using the definition ( p.l9|) of V{r). It is 
conveniently partitioned as 

H = Hi + H2 + Hi, (3.6) 

with 

H, = J^-^, (3.7) 

= - , (3.8) 
2m 47rr2 

^^ = :r^ (3-9) 

47r|ri — 47rr2 

We note here that 0j(r2) is the ground state wave function for the Hamiltonian H2 defined 
in ( p.8|) . This is the reason why we split the total Hamiltonian in the manner shown above. 
We shall find this choice makes the evaluation of the direct term easier. When we consider 
the exchange term in appendix D, the definitions of Hi and H2 are switched since there ri 
and r2 are interchanged. 



A. Direct Terms; Comparison Potential Choice 



We now examine the direct term T^. Expanding Eq. 
leading order term for Td, 



in powers of Hj, we get the 



2e' 



47r|ri— r2| Anri 



-viri) 



1 



_ Hi + H2 - E - Ef - le 



|ri=0, i) 



2e' 



-v{ri) 



Hi+Ei-E-Ef-it 



|ri=0), (3.10) 



47r|ri— r2| Auri 

where we have used H2\i) = E\i). Since the final ion S-states contribute the major part of 
the total decay rate, we shall first consider only the cases where the final ionic state is in 
S-state. We shall choose the comparison potential to be given by 

,2 . 1 2e2 



v{r) 



e 
47r 



J r — r 



Airr 



(3.11; 
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where0 



(3.12) 



It is easy to see that with this choice the leading order term ( p.lOp vanishes. 

With the choice ( p.ll| ) of the comparison potential, we show in Appendix A that, including 
terms up to order 77^, 



.(0) 



-e^v r(i _ 2tr]')e 



(3.13) 



which is a modification of the result in reference The phase G, which is irrelevant in 
the beta decay problem, is defined in Appendix A. The shifted momentum p' is related to 



fi 



{dh)pfi 



a 



Airr ao 



Vfi 



(3.14) 



via 



P 



— 2mVfi = Pyl — 2rfvf^ 



(3.15) 



and 



rj' = am / p . 



(3.16) 



Note that Vfi is the potential energy at the origin produced by the charge distribution (|3.12|) , 
which involves both the initial and final states, while Vfi is a. dimensionless potential. This 
result differs from the correction derived in the literature which uses the charge density 
of the initial state, not our "transition density" pji. 



^'^This density is always real because the wave functions for the initial and final states can be written 
as a real function multiplied by a constant phase factor which is canceled by the normalization factor 
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B. Exchange Terms 



The third term T3 in Eq. ( |2.27| ) and the term Tg in Eq. (3^) are both exchange terms. 



They are examined in Appendix D, where the resuh ( pi8| ) shows that the exchange am- 
phtudes are of order rj^. This contradicts a result appearing in the hterature |14], where 
the leading exchange amplitudes are claimed to be of order t]^. Though the exchange am- 
plitudes are of order rj^, they only contribute to the decay rates at order 77^ since they are 
relatively imaginary to the leading sudden approximation amplitude. The details are given 
in Appendix D. 

C. Order t]^ Corrections 

The first correction to this leading order direct term Td is obtained by keeping one more 



term while expanding Eq. (3^) in powers of H 



If- 



n = {P,f\ 



^2 2e2 



47r|ri — Atcti 



- f (n) 



This leading term of is of order 77^. Therefore, to the second order of 77, we may do 
following approximations. First, it is valid to drop Ef and H2 in the denominators relative 
to energy E since they are of ordei0 rj'^E. Secondly, we can approximate the beta ray wave 



^^We provide here more justification for dropping the Hamiltonian H2. Since the wave function 
(/)j(r2) is the ground state eigenfunction of H2, one can replace H2 in the second denominator by 
Ei. For the first denominator, we can imagine inserting a complete set of eigenstates of H2 just 
before Hj in the second line of Eq. ( 3.17| ). Since the wave functions (/>t(r2) and 0j(r2) are slowly 



varying functions, even after multiplied by the Coulomb interaction factor l/|ri — r2|, they have 
little overlap with the eigenfunctions of H2 with energies much higher than the atomic energy due 
to the highly oscillating feature of these high energy eigenfunctions. Therefore, in the sum over 
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function 0p(ri) by the plane wave exp(— ip-ri) or equivalently treat (p| as a free particle state 
with momentum p since their difference is of order t]. Finally, the Hamiltonian of electron 1 
Hi may be replaced by its free Hamiltonian Hq. This is justified because Hi describes the 
motion of the beta ray, thus the Coulomb interaction provides only higher order corrections 
in the rj parameter to its propagation. With these approximations, to order r]'^, we get on 
using Eq. (gj), 

e2 2e2 



T]^ /(rf^r2)0Kr2)0.(r2)(p 



Ho-E-ie 



47r|ri - 
/ 1 
An \ Iri - 



Anri 

--) 



v{ri) 
1 



Ho-E-ie 



|0) 



(3.18) 



In Appendix B, we find that J/j is given by 



Jfi = iffi - 1 - , 



(3.19) 



with 



r(p) = Tf,^^e^^T{l - 2tr])e- 



i0 



1-^(1- ff^ + Kj. 



(3.20) 



Combining these second order corrections ( p.20|) with value of the wave function 0p(r) at 
the origin ( p.l3| ) gives the the atomic matrix element T(p) to order rf: 



(3.21) 



Taking the square of this matrix element, using the definition ( |3.15| ) of p', dropping higher 
order terms in 77, and recalling the differential decay rate formula ( |2.36| ) yields the exclusive 
differential decay rate to an ionic S-state: 



the eigenstates of H2, the main contribution comes from summing over the excitations with energy 
being of the order of the atomic energy. Thus, in both denominators, H2 is of the order of the 
atomic energy. 
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^ = ^F(2,E)p|T,r(Q,-Ef |T,„p [l-2r^^(l + i^4 , (3-22) 

where F{2,E) is the Fermi function defined in Eq. ( p..?!) . Note that the term in Eq. (|3.21| ) 
involving Vfi cancels the term from expanding the momentum ratio p'/p. 

Using the wave functions listed at the end of Appendix B, it is easy to compute 

, ,5 fl6V2\^ 512 

Tij2 = — — = — ~ 0.7023 , 3.23) 

I 1^1 \^ 27 / 729 ^ ^ 

\T2i\^= f-^) =7 = 0.25, (3.24) 



and 

|T3.P=(-^^1 =^^ 0.0127, (3.25) 



144^6 \^ 2^3^ 



which demonstrates the dominance of the first two exclusive decay rates. It is also straight 
forward to compute Kfi for a transition to any particular final state. For example, we show 
in Appendix B that for the transitions to the first three low energy S-states, 

975, (3.26) 
^2. = - -J ^ -0.3025 , (3.27) 
K3.= f--^U -0.7902. (3.28) 
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12 


-Ih 




1045^ 


12 


648 J 



IV. SUMMING OVER ATOMIC S-WAVE STATES 



We shall denote by Fq the contribution to the inclusive decay rate from the cases where 
the final ion is in an S-wave. According to Eq. (|3.22| ), this is given by 



dE 27r3 



F(2,i?)p|T^P {Qf-Ef\Tf,\'' l-2rf{l + kf.. 



(4.1) 



/</n 



where the upper limit /max of the summation corresponds to the final ionic states with energy 
Ef being -Emax = Q + Ei — E since the total energy must be conserved. In view of Eq. (|2.37|) , 



23 



{Qf - Ef = {Q-Ef + 2{Q- E){Ei - Ef) + {E^ - Ej 



(4.2) 



we need to do sums weighted by {Ei — Ef)"', with n = 0, 1, 2 . The result will be expressed 
as the sudden approximation spectrum plus the r]'^ order correction: 



dTn m 



dE 2tt^ 



E{2,E)p\Tp 



'l-2rl')P{Q-E)-\c{Q-E) 



(4.3) 



We shall first examine the sums corresponding to the sudden approximation which defines 
P{Q — E) and then the sums involving Kfi which defines C{Q—E). 

A. Sudden Approximation Terms 



To do the sum over final states which do not include energies higher than 

_ 

7, = -^max = Q + Ei — E , 

2m 



we write 



E l/)(/l = i- / 



k)(k| 



Using the squared matrix element calculated in Appendix C, Eq. 

2567r2 1 76 



-87 cot ■'7 



(4.4) 



(4.5) 



(4.6) 



1 - fc3 (1 + ^2)4 

where 7 = 1/kao, we can evaluate the needed sums appearing in the sudden approxima- 
as 



tion (g]T 



f~r Jk>K (ZTT)'^ 

/<./max ^ ' 



k>K (27r) 
32 /-V^ao 



(^\{H,-Hjn)-- dj 



l_g-47r7 (1+72)4 

where E^ is the energy of the ionized electron, 



r 



(1+7T , (4.7) 



Ek = — 
2m 



(4i 
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and % is the Rydberg constant, 



1 



13.6 eV. 



(4.9) 



The first term in Eq. (|4.7|) is the closure approximation result, and it can be calculated easily 
by using the ground state wave function of hydrogen atom: 



1, n = 0, 
2Ry, n=l, 
8 %^ n = 2. 



(4.10) 



The sudden approximation spectrum is proportional to the sum 

P{Q-E)^ Y: \TM'{Qf-E)\ 

/</m ax 



(4.ii; 



which may be expressed in terms of the closure approximation result plus the correction 
to the closure part due to the fact that the summation does not include the final states 
with energy higher than -Emax = K'^/2m. Since K is related to Q — by the relation 
-E'max = K'^ /2m = Q — E + Ei, we write P as a function of variable Q — E for the convenience 
of later usage. Explicitly, using the closure results ( [4.10| ) above, we have 



P{Q - E) = {Q - E + 2%)2 + 4V + S{Q-E) 



(4.12) 



where S{Q—E) denotes the correction to this closure approximation result. In view of 
Eq. ( |4.7| ), this correction comes from the 7 integrals, which can be evaluated numerically. 

The function Pi{E — Q) defined by Eq. ( |1.14| ) in Section I for the purpose of estimating 
the atomic effects in neutrino mass measurements may be obtained immediately by setting 
n = in Eq. (|3) : 



32 rl/Kao 4ti^ 

P^{Q-E) = l-- dj- ^7T— ^ 

n Jo 1 — e ^''^ (1 + Y) 



^ g— 87 cot -'-7 



(4.13) 



This expression enables a numerical evaluation of the function Pi{Q — E). 

Before displaying the numerical results, we consider the asymptotic behavior of S{Q~E) 
and Pi{Q — E) as Ka^ becomes large. For Kqq ^ 1, we study the asymptotic expansion 
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of the 7 integrals in Eq. ( |4.7| ). This may be accomphshed by expanding the integrands in 
powers of 7 since the integral takes the value of the integrand in the interval (0, 1/ Kao). One 
can then do these integrals easily. Keeping only the three leading terms in the expansion, 
the large Ka^ asymptotic form for S{Q—E) is 



S{Q-E) ^ ^ 



1 5 /327r 32 \ 1 



By\ (4.14) 



nKao {Ka^f V 21 It: J {Kaof _ 
As one of the intermediate steps above, the leading large Ka^ behavior of Pi {Q — E) is 

32 1 

We show in Fig. 1 the numerical result of S{Q—E) as a function of -Emax = K'^/2m = 
Q + Ei — E. The asymptotic form ( [4.14| ) is also shown in Fig. 1, where one finds that 



it describes S{Q—E) with a good accuracy until KaQ is less than 5, which is expected by 
observing that the coefficients of the three terms in the asymptotic form ( ^.14| ) are all of 
order 1. The numerical result of function Pi{Q — E) is shown in Fig. 2 together with its 
asymptotic form ( ^.15| ) as functions of -Emax- 



B. Kfi Corrections 

In this section, we shall deal with the contribution to the inclusive differential decay rate 
involving Kfi. Recalling the differential decay rate formula ( [1.1| ), the sums which we need 



to consider involve the combination \Tfi\'^Kfi. Therefore, it is convenient to define 

K{e)^{z\r-V){f\r^\t), (4.16) 



which is related to Kfi via 



\Tf.\'Kf. = \K"{0), (4.17) 



in view of Eq. (|3.20|) . Correspondingly, the sum 



Kn{e)^ E {^\r-V){f\r'\^){E.-Efr, (4.18) 

/</max 
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FIG. 1.: Curves for S{Q—E) — the correction to the closure part of the sudden approxi- 
mation — as a function of the energy threshold of the final ionic state, -Emax = K"^ /2m = 
Q + Ei — E and its asymptotic behavior for large Q—E. The solid curve is the plot of numer- 
ical result of S{Q—E), the lower and dashed curve is the high energy asymptotic behavior 
given in Eq. ( [4.14|) . The units are for the vertical axis and for the horizontal axis. 
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FIG. 2.: Curves for Pi{Q — E) — a function measuring the change of the spectrum due to 
per unit change of the neutrino mass squared — as a function of the energy threshold of the 
final ionic state, -Emax = K'^/2m = Q + Ei — E and its asymptotic behavior for large Q—E. 
The solid curve is the plot of numerical result of Pi{Q — E), the lower and dashed curve is 
the high energy asymptotic behavior given in Eq. ( glSl) . The units are 1 for the vertical axis 
and Ify for the horizontal axis. 
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can produce the desired sums through 



E \Tf,\'KfaE,-Efr = -K{0). 

f<fm ax 



(4.19) 



Since the spectrum is corrected by the combination I]/</max I'^M'^^fiiQf ~ i ™- view of 
the relation ( [4.19| ), the contribution to the spectrum due to Kfi can be expressed in terms 
of the quantity C{Q—E) defined by 



C{Q-E) = {Q- EYK'^iO) + 2{Q - E)K'l{0) + K'^{0) 



(4.20) 



The goal of this section is to show the main steps of calculating C{Q—E) and display the 
numerical results of C{Q—E). 

For n = 0, the completeness relation 

(rf^k) 



and the e derivatives give 



Km 



^o(e) 



k>K (27r)3 



1 



(4.21) 



l/{Kao) 



d-i-i-^K"{<S). 



(4.22) 



lk>K (27r)3" 27r2a|] 
With this, Kq may be then directly evaluated by exploiting the numerical results for K"{{]) 
obtained in Appendix C. 

For = 1, we may write similarly 



K,{e) = {^\{H,r^-r^Hf)r-^\{)- j ^^,K{e){E,-E, 

Jk>K [ZTiy 



(4.23) 



The first, closure-approximation term is readily evaluated: 



{i\{Hir' -r'Hf)r-'\i) = {i\ 



ie / r r 

47rr 2mr2 2m \ ^ 



\^) 



2(l + e2)Ry. 



(4.24) 



Therefore, 



4 + 
4 + 



(rf^k) 



k>K (27r)3 

1 r'^HKao) 



K"mi+i'') 



Jo 



rf77-6(l + 72)K"(0) 



(4.25) 
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which enables again a simple numerical calculation of K"{0) upon using the result of K"{0). 
For the case n = 2, since the closure approximation result diverges, we shall use a 



different formalism. It is convenient to separate the sum in the definition ( |4.18| ), and thus 
the definition for (0), into two parts: the contribution from summing the final ionic bound 
states denoted by 

K,^ = Y.K"mE,-Ef)\ (4.26) 

b 

and the contribution from the continuum with energy under the threshold K"^ /2m denoted 
by 

Jk<K [iT^y^ I'K'^ Jl/Kao 

The sum of and produces (0)- The second part can again be readily evaluated 
by using the numerical results for K"{Q) given in Appendix C, and this Appendix also 
provides the evaluation 

i^zb = 32.26 V- (4.28) 

The curve for C{Q—E) as a function of the threshold energy -Emax is shown in Fig. 3. 
We find that a quadratic function Cfi{Q—E), 

Cfi (Q-E) = 33.6 V + 4.92 ^r E^^, + 0.0148 

= 28.7 V + 4.89 Tfy{Q-E) + 0.0148 (Q - Ef , (4.29) 

describes C{Q—E) for -Emax in the range — 64% (0 — 870 eV) with a good accuracy. To 
analyze how the atomic effects change the experimental data fitting, we fit C{Q—E) to a 
linear combination of the functions P{Q—E), P'{Q—E), and Pi{Q—E) in the same range of 
the spectrum and get 

C{Q-E) ^Cf2{Q-E) =0.0U8P{Q-E) + 2A2ByP'{Q-E) + 18.9By^Pi{Q-E) . (4.30) 

We show C{Q—E) and its two fitting formula Cfi{Q—E) and Cf2{Q—E) simultaneously in 
Fig. 4. The three curves can barely be distinguished. The errors due to the fitting are shown 
in Fig. 5 for Cfi{Q—E) and Fig. 6 for Cf2{Q—E) respectively. 
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FIG. 3.: Curve for C{Q—E) — the correction to the spectrum due to Kfi as a function of 
the threshold energy of the final ionic state E^^x — K'^/2m — Q + Ei — E; the units are By^ 
for the vertical axis and By for the horizontal axis. 
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Emax 



FIG. 4.: Curves for C{Q—E) and its two fitting forms Cfi{Q—E) and Cf2{Q—E) which can 
barely be distinguished. The units are the same as in previous figure, By^ for the vertical 
axis and By for the horizontal axis. 
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FIG. 5.: Curve for the difference Cf\{Q—E) — C{Q—E) as a function of ii^max = Q + Ei — E. 
The units are the same as in previous figure, for the vertical axis and By for the horizontal 
axis. 



33 




10 20 30 40 50 60 



E max 



FIG. 6.: Curve for the difference Cf2{Q—E) — C{Q—E) as a function of ii^max = Q + Ei — E. 
The units are the same as in previous figure, for the vertical axis and By for the horizontal 
axis. 
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V. FINAL ATOMIC STATES WITH NON-ZERO ANGULAR MOMENTUM 



So far, we have considered only the final ion state (/| being in an S-state. To obtain 
the inclusive decay rate, we also need to consider the case where the final ion states are 
not S-states. For these cases, the first term Ti, the amplitude of the sudden approximation, 
vanishes. Therefore, to order iff, as far as the decay rate is concerned, it is sufficient to 
calculate the amplitude to order r/, which comes from the leading term of T^. Recalling 
expression for this leading term T^, we have, to order 77, 



T°^-(p||(rf3r2)0}(r2)0.(r2) 



-v[ri) 



^ 'ri=0). (5.1) 



Hi+Ei-E-Ef-ie 



47r|ri— Attti 

The last two terms in the square brackets vanish upon integrating over the solid angle of r2 
since now </>/(r2) contains only higher partial waves with Z > 1. Hence 

T° ^ - / {d'r,)<j>}{r,)<j>,{r,){p\—^ r ^ ■ |0) , (5.2) 

J ■' 47r|ri — Hq — E — le 

where we have neglected Ef compared with E and replaced Hi by the free Hamiltonian 
Hq. Suppose that final state (/| has the angular momentum quantum number {l,m), or 
equivalently 

0}(r)=i?},(r)y^:„(r), (5.3) 
with Rfi{r) the radial wave function. In Appendix B, we find that, to the leading order, 

Ayr fOO 

^ */(7^^^'-(P) I dr2rlR},{r2)Ur2) • (5.4) 

The differential decay rate involves the angular average of the square of the spherical 
harmonic function which appears in Eq. (|5.4|) , 

and thus, in view of Eq. ( p.36| ), the differential decay rate to states with specific energy and 
angular momentum I is given by 
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dr 



dr2r^R)i{r2)<i)i{r-2) 



(5.6) 



dE 27r3^' ' ' /2(/ + i)2 

where we have done the summation over the magnetic quantum number m which generates 
the factor 2/ + 1. 

Summing over I and / produces the differential decay rate for the final ^He^ ion having 
nonzero angular momentum. Using the expansion ( |4.2| ) we encounter the sumslll 

47r(2/ + 1) 



roo ^ 



(5.7) 



with n = 0, 1, 2. To facilitate the calculation, we define 



M{k) = {E,-EkfJ2 



47r(2/ + 1) 



drr^Rli{r)(J)i(r] 



(5. 



where Rki{r) is the radial wave function of the final unbound ion state with energy Ek = 
k"^ /{2m) and angular momentum /. M{k) is calculated in Appendix C. We can consequently 
calculate Mn for n = 0, 1, 2. 

For n = 0, exploiting the completeness relation 



Y.R),{r')Rfi{r) = -^5{r'-r) 
f ^ 



(5.9) 



we can write 



V Pil+iy Jo 

roo 2A;2 



drr'^Rki(r)(f)i(r] 



^~ jj^ dk—M{k){Ei-Ek 



(5.10) 



where the sum 



2/ + 1 



E 



(5.ii; 



and the unit norm of the initial wave function 0i(r) have been used. This expression enables 
a numerical evaluation of Mq. 



Here the sum over / should still be understood as the sum with the upper bound /max as before. 
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To calculate Mi and M2, we write M„ as the sum of two parts — M„b the part coming 
from summing over the final bound ionic states and Mnc the part coming from summing over 
the final continuous ionic states, 

M„ = M„b + M„c. (5.12) 

Following the definition of Mnc, we have 

Mnc= dk M{k){E.,-EkT-\ (5.13) 

which may be easily calculated numerically. For M„b, we shall not investigate it numerically; 
instead, we estimate it to a good accuracy. To facilitate the notation, we define L{Q—E) by 

L{Q-E) =Mo{Q- Ef + 2Mie {Q ~ E) + M^c , (5.14) 

which enables us to write the differential inclusive decay rate for the final ion state having 
nonzero angular momentum as 

^^T.^ = ^sP\Tp\'v' [L{Q-E)+2MM-E)+M,,] . (5.15) 

Upon using the numerical result of M{k) obtained in Appendix C, one may obtain Mq, Mic, 
and M2C by numerically performing the integrals over k in expressions (|5.10|) and (|5.13|) . 
Definition (|5.14|) thus gives a numerical result for L{Q—E) which we display as a function 
of -Emax = K'^/{2m) in Fig. 7. We find also a quadratic function Lfi{Q—E) 

L/i(Q-E) = 5.90%' - 1.03 Ry^max + 0.947 

= 7.88 V - 2.93 Ry (Q-E)- 0.947 {Q - Ef (5.16) 

which fits L{Q—E) for -Emax in the range — 64 (0 — 870 eV) with a fairly good accuracy. 
Like in previous section, we fit L{Q—E) to a linear combination of functions P{Q—E), 
P'{Q—E), and Pi{Q—E) to generate the fitting formula 

L{Q-E) ^ Lf2{Q~E) = 0M8P{Q-E) - 3.39ByP'{Q-E) + 15.7By'^Pi{Q-E) . (5.17) 
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FIG. 7.: Numerical curve for L{Q—E) — the correction to the spectrum due to atomic effect 
for non-S-wave final ionic states — as a function of E'max — Q + Ei — E; vertical axis has the 
unit By^ and horizontal axis has the unit By. 
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FIG. 8.: Numerical curve for L{Q—E) and its two fitting formula Lfi{Q—E) and Lj^iQ—E) 
as functions of E'max — Q + Ei — E; vertical axis has the unit By^a^ and horizontal axis has 
the unit By. 
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FIG. 9.: Numerical curve for the difference between Lfi{Q—E) — L[Q—E) as a function of 
-E'max — Q + Ei — E, vertical axis has the unit By^ and horizontal axis has the unit By. 
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FIG. 10.: Numerical curve for the difference between Lf2{Q—E) — L[Q—E) as a function of 
Eraax — Q + Ei — E, Vertical axis has the unit By^ and horizontal axis has the unit By. 
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We show L{Q—E) and its two fitting functions Lfi{Q—E) and Lf2{Q—E) together in Fig. 8. 
They can hardly be distinguished. The errors caused by the fitting are shown in Fig. 9 for 
Lfi{Q—E) and Fig. 10 for Lf2{Q—E) respectively. 

We now turn to provide numerical bounds for Mib and M2h- The energy of an ionic state 
with principal quantum number Uf is given by 



Since the lowest energy state with nonzero angular momentum {uf — 2) has the energy Ei, it 
does not contribute to the summation defining M„b for n = 1, 2 due to the factor (E^ — EfY. 
For the other states, since Uf > 2, Ef > 4£^j/9, and hence 



where Mq^^ is defined by excluding the lowest energy state {nf—2, l—l) in the {rifj) summa- 
tion of the definition of Mob, ^-c. 




(5.18) 




(5.19) 



Therefore, according the definition of M^h, we have 




(5.20) 




(5.21) 



Evaluating the matrix element. 




(5.22) 



and using the numerical result 



Mob = 0.659 , 



(5.23) 



gives the explicit bounds: 



-0.05 Ify > Mib > -0.1%, 



0.03%^ < M2b < 0.1 IV^ 



(5.24) 
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VI. INCLUSIVE DECAY RATE 



Adding the differential decay rates (O) and ( 5.15| ), we find that the inclusive differential 
decay rate is given by 



dE dE ^ dE 
m 



F(2, E) p \Tf,\'{l-2r]') P{Q - E) + r^'R{Q-E)] 



27r3 

where we have defined the correction term R{Q—E) by 

C{Q-E) 



(6.1; 



RiQ-E) 



+ L{Q-E) + 2Mib(g -E) + M2b : 



(6.2) 



and dropped higher order terms. Here C{Q—E) and L{Q—E) are evaluated numerically; Mib 
and are very small contributions bounded by Eq. ( ^.24|) . The term containing R{Q—E) 
in the second line of Eq. (|6.1| ) represents the correction to the sudden approximation result 
due to the atomic effect. 



A. Comparing with previous results 



We now compare our result ( |6.1| ) with previous results |rT] , P^ . To reproduce the previous 
results, two approximations must be made, as mentioned in Section I. Under the "uniform 
phase space factor approximation" , we have 



C{Q-E) ~ (Q - ^)X(O) 
L{Q-E) ~ {Q-E^Mo, 



(6.3) 



where E is some average beta ray energy which differs from the beta ray energy E by an 
amount of order I^Jll Under the closure approximation, Eq. (f4.22|) and Eq. (|5.10|) read 



^^For example, in reference E is chosen to be — 3Ify. The details of the choice of E does 
not matter, since there the order rj'^By/(Q — E) terms are omitted. 
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Kq{0) ~ and Mq ~ 1 respectively. With these approximations, the definition (|6.2|) gives 
R{Q—E) ~ (Q — E)"^, which imphes that the result ( |6.1| ) above reduces to a modified sudden 
approximation result: 

^ = £^F(2,Ey|T,pP(g-E), (6.4) 

where terms of order 'rfS{Q—E) and "rf^"^ have been discarded with the approximations 
we are considering. The result ( |6.4|) agrees with previous results [irT||l2| , with the modified 
momentum p' being defined by 

The momentum p' equals the momentum of an emitted beta electron whose energy at short 
distances is modified by the repulsive Coulomb interaction energy with the original electron 
bound in the tritium atom. This modification which changes p to p' was found by Rose a 
long time ago |T^. 

B. Estimation of the atomic effects in the neutrino mass determination 

To estimate how the atomic effects change the neutrino mass squared parameter, we first 
recall that the spectrum used in the experimental data analysis is the sudden approximation 
result 



AF{2,E)p Y: Wf,{Qf-E)^{Qf-EY-ml. (6.6) 

""-^ /</m ax 

The terms in the right hand side of Eq. (|6.6| ) require some explanation. The first factor 
F(2, E) is the usual Fermi functional with the nucleus charge Z = 2, and p is the momentum 
of the beta ray. Wji is the transition probability for the initial tritium state decaying to 
the state / of the final '^He^ ion. In real experiments, molecular tritium is used. Therefore, 



14 



Usually a relativistic Fermi function is used to take care of the dominant relativistic correction. 
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Wfi is the square of the matrix element for the molecular state. Since only atomic tritium 
is considered in this article, we have made the replacement Wfi = Finally, is the 

neutrino mass squared. Though one focuses on measuring the beta ray spectrum near the 
end point to probe the neutrino mass, most of the data obtained in the experiments are in 
the range where {E — Q)^ ^ since the decay rate is tiny at the very end of the spectrum. 
Expanding the square root, we get 



exp 



dE 



AF{2,E)p Y: \TM'[{Qf-Ef--ml 



A F{2,E) p 



P{Q-E) - -mlP,{Q-E) 



(6.7) 



where we have used the definitions ( |1 . 1 1|) and ( p..l4| ) for the functions P{Q—E) and Pi{Q—E) 
which are evaluated numerically in section IV. 

In real experimental data analysis, the parameters A, Q, and ml are determined by 
comparing the spectrum (|6.7|) with the measured spectrum. Therefore, any theoretical cor- 
rection to the spectrum ( |6.7|) has the effect of changing the parameters A A = A + AA, 
Q ^ Q = Q -\- AQ, and ml Ami thai the theoretical correction to the spectrum may 



be included by using an effective spectrum described by the same form (|6.7|) but with the 
effective parameters A, Q, and Ami. The change of the spectrum ( |6.7| ) 



\ dE 



A F{2,E) p 



A A 1 

— P{Q-E) + AQP'{Q-E) - -AmlP,{Q-E) 



{61 



mimics the corresponding theoretical correction to the spectrum. 

We now estimate AQ, Ami corresponding to the correction due to the atomic effect 
which accounts for the interaction between the beta ray and the electron of the ^He'*' ion. 
This we shall do by requiring the correction R{Q—E) in Eq. ( |6.1| ) be mimicked by a linear 
combination of P{Q—E), P'{Q—E), and Pi{Q—E) as appearing in the right hand side of 
Eq. ( |6.8| ). We can neglect the small parameters Mib and since they are bound by 
Eq. ( |5.24| ). Since the region important for the neutrino mass measurement goes from the the 
beta ray end point to approximately 59 — 74Ry (800 — 1000 eV) below end point 
the previous fitting formula Cf2{Q—E) and Lf2{Q—E) for the energy -Emax range — 64Ry 



45 



may be used here. We shall discuss the sensitivity to the range of the energy of the result 
later. Replacing C{Q—E) and L{Q—E) with their fitting formulas ( [4.3CI| ) and ( p. IT)) , the 
definition (16.21) reads 



R{Q-E) ^ OMP{Q-E) - A.QBy P'{Q-E) + Q.3By^ Pi{Q-E) . (6.9) 
Inserting this fitting form of R{Q—E) into Eq. ( |6.1| ) yields 

^^^,F{2,E)p\mi-imv-) 

X 'P{Q-E) - 4.6 T]^ By P'{Q-E) + 6.3 r/^ V Pi{Q-E)] , (6.10) 

where we have shifted the argument of P to absorb the second term which causes an negligible 
error of order 0{ri'^). In view of the argument above, specifically Eq. ( p.8|) , Eq. ( |6.10|) shows 
that the atomic effect gives a correction to neutrino mass squared of 

Ami ^ -12.6 ri^By^ -1.7 eV^ , (6-11) 

and the endpoint changes by 

AQ ^ -4.677% ^ -0.047 eV. (6.12) 

We now examine the sensitivity of the result to the energy range i?max used in the fit. Since 
AQ is tiny, we shall only consider how Ami depends on the range in which we do the fit. 
We fit R{Q—E) in various energy -Emax ranges and display the corresponding ml in Table 1 
shown in the introduction. The neutrino mass squared has basically a linear dependence on 
the energy range where we do the fit. Increasing the energy range by each 5 Ry causes Ami 
to decrease by — 0.4eV^. The atomic effect changes the neutrino mass squared parameter on 
the order of a few eV^. It is not a big effect. 



VII. CONCLUSION 



We have developed a systematic expansion for the tritium beta decay amplitude in the 
Coulomb parameter rj. By choosing a convenient comparison potential, one can avoid the 
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infrared divergences due to the long range characteristic of the Coulomb force. Both the 
exclusive and the inclusive decaying rates are calculated to order rf. The inclusive decay 
rate agrees with previous results. The estimation on how this order rf correction affects the 
neutrino mass squared parameter is provided. We find that the effect is small and does not 
suffice to explain the mysterious negative electron anti-neutrino mass squared obtained in 
modern experimental data analysis ||l],H^^J^,^] . We also remark that the order iff correction 
to the spectrum is small and can not provide any explanation for the anomalous structure 
in the beta decay spectrum in the last 55 eV closest to the end point presented in the last 
article of reference 0. 
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APPENDIX A: WAVE FUNCTION CALCULATION 



We shall find the value of the wave function 0p(r) at the origin up to and including terms 
of order rf. Since 0p(O) involves only the S-wave component, we may replace 

0p(r) ^ —Up{r) , (Al) 
pr 

where Up{r) obeys the S-wave radial Schrodinger equation 



+ v{r) - E ) Up{r) = . (A2) 



[ 2m dr"^ 

The radial wave function Up{r) vanishes at the origin and obeys the asymptotic boundary 
condition 



^ ^ ^i'pr+ir]ln(2pr) ^—ipr~iT]ln(2pr)—2i5o 



(A3) 



where 5o is the S-wave phase shift. This boundary condition contains a non-trivial phase 
structure because, at large distances, the particle moves in the long-range Coulomb field of 
unit charge. 

It is convenient to divide space into two regions: < r < and < r < oo, with 
being an intermediate distance between the Bohr radius Oq and the de Broglie wave length 
= TjaQ. We take to be of the order 

(r,/ao)' ~ O(r/i+0 , (A4) 

with ^ < C < 1 [e. g. Tc = rf^'^ao). As we shall see, there are two reasons for this choice. One 
reason is that in the region where < r < rc, we have r <^ ao which enables us to expand 
the expression ( |3.11[ ) for v{r) as 

v(r) = h Vfi + Oie^rfp^r'^) , (A5) 

A-Kr 

where 

2e2 



Vfi = lim 



f (r) + 



47rr 



(^'r)p^,,(r)^. (A6) 
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The expansion (|A5|) of Av{r) does not contain a term linear in r. Such a term would 
correspond to a charge distribution — V^r ~ 1/r which has an unphysical singularity at 
the origin. The correction to the first two terms in the expansion ( [A^ ) comes from the 
consideration that, for any reasonable comparison potential, the characteristic length scale 
for Av{r) to vary is the Bohr radius a^, and so the leading correction is of the order 



/ r \ ^ 
r KaoJ 



(A7) 



Inserting the expansion (|A5|) of the potential into the radial Schrodinger equation ([A^), we 
see that the "interior" solution Up{r) in the region < r < is a constant C times the 
Coulomb S-state radial wave function with charge 2 and energy Ep/ = E — Vfi. This radial 
wave function involves the shifted momentum 



p' = Jp^ — 2mVfi 



(A8) 



and the correspondingly altered Coulomb parameter 



2ri' = 2am/p' . 



(A9) 



In the region r < r^, the higher order terms in the expansion ( |A5D give rise to corrections, 
compared with the energy E of the beta ray, of order 



2 2 
6 'T'j 1" 



E 



(AlO) 



for our choice ( |A^) of r^. Therefore, this Coulomb wave function with the shifted momentum 
p' obeys the Schrodinger equation ( |A2|) with an error which is less than order r]^. As is well 
known ||16[, this wave function gives the limits. 



r : Up{r) CpVe^'''r(l + 2ir]') 



(All) 



and 



r — > oo : Up{r) ~ C 



At]'' 
pr 



-1/4 



1 ip'r+2i'r}' In p'r _i_ 

2t\ 



(A12) 
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where the eUipsis ■ ■ ■ stands for the incoming wave contribution which involves exp{—i{p'r + 
2r)'ln2p'r)}. Note that, for r ~ r^, the asymptotic expansion ([A12[) is vahd because pr ^ 1 



in this region. This is the other reason for the choice (^1) of r, 



In the region where r > Tc, we can get the asymptotic form of the wave function Up{r) by 
iterating Eq. ( [X^ ) starting with the hmiting behavior ( |A3D of Up{r) for r — oo. Matching this 
"exterior" solution with the previous "interior" solution at r ~ rc determines the constant 
C. The two linearly independent solutions at large radius r correspond to outgoing and 
incoming waves, 

Up{r) = (r) + (r) . (A13) 

Since we have only one constant C to be determined, it is sufficient to match the outgoing 
wave part u^\r) to connect the solutions in the two regions. We write this part of the 
exterior solution as the W.K.B. approximate solution times an arbitrary function. 



1 



Here S{r) is defined by 

dS{r) 
dr 

with the boundary condition that 



E 



E — vir) 



1/4 



;iS(r) (A14) 



2m[E -v{r-)], (A15) 



S'(r) —> pr + 7] In 2pr , asr-^cxo, (A16) 

which is consistent with the long-range limit of v{r) in Eq. ( |A15| ). Thus, the boundary con- 
dition (^) is obeyed by requiring that w':^\r) — > 1 as r — > oo. Iterating the solution ( |A14| ) 



in the Schrodinger equation (^) yields the asymptotic expansion of w^'^^(r). For example, 
the first iteration gives 

/ , N , , im dv(r) 



When r ~ Tc, which is in the asymptotic region, the form ( [A14| ) for u'^\r) with w^+^(r) = 1 



is appropriate for our discussion. In the region r ~ r^, the first factor containing the power 
1/4 in the solution ( |A14|) , by using the expansion (|A5|) , can be written as 
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which matches the corresponding factor in the asymptotic expansion ( |A12| ) of the interior 
solution. The phase S{r), in this region, to order t]^, is given by 

S{r) ^ p'r + 2r]' In 2p'r + Q , (A19) 

where is a constant phase. The fact that is a constant may be justified by observing 
that 



[S{r) — p'r — 2r( In 2p'r\ = \j2m{E — v{r)) 



-p - 



dr * 2p^ r 

p 2p^ r r, ' 



where we have used the definition of S{r), expanded the square root in powers of v{r), and 
replaced v{r) by its expansion (|A5|) .E1 Therefore, in the region where r ~ Vc, the leading 
terms in the solutions ( [A14| ) and ( |A12|) match to the order t]^. We did not include, in 



the asymptotic expansions ([A12| ) and ( |A14|) , sub-leading terms in 1/r such as ///(pr)^ and 
ri/{pr)^ which are of the order 0{7f) and 0{if') in the matching region. These terms must 
also match to the order 7f because we are solving the same differential equation by using 
different approaches with both having an accuracy of order For example, the first 

correction in Eq. ( [A 171 ) shows the matching of terms of order r]/ (pr)"^, if we recall the leading 



behavior of v{r) at the region r ~ and include the term of order ri/{prY in the asymptotic 
behavior (|A12D . 



Requiring that the outgoing wave in the asymptotic form ( A12 ) to be the same as that 
of M^+^(r) in Eq. (|A14|) thus gives 



15 



The phase may be exphcitly determined to 0{rf) by integrating Eq. (A15) inward from the 



hmit ( |A16D . 

^^This remark shows that we have presented more detail of how the matching works than is actuahy 
necessary. This we did for the sake of clarity. 
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V p' 



(A21) 



Consequently, inserting this expression for C into Eq. ( [All| ) and using Eq. ( |70] ) gives, in- 
cluding terms up to order rj^, 



()p(0) =e"'''r(l + 2zV)W-e'® 

VP 



(A22) 



The phase G is, of course, irrelevant since only the absolute value of the reduced matrix 
element occurs in the decay rate. 

APPENDIX B: SECOND ORDER CORRECTIONS 



1. Evaluation of J 



fi 



We now calculate Jfi defined by Eq. ( p.l^ ). Defining 



^p(r) = (p| 



1 1 



|ri— r| ri 
with Av{r) being given by 



-Av(ri) 



1 1 



HQ—E—ie \|ri— r| r I HQ—E—ie 



-|ri=0) (Bl) 



expresses Jfi as 



The representation 



—Av{r) = v{r) + - — , 
47r Airr 



dte-'tiHo-E) 



Hq — E — it Jo 



expresses 



roc rco 

Lp(r) = - / dt dt'{p 
Jo Jo 



1 1 
Iri— r| ri 



-At>(ri) e 



-it{Ho-E) 



_j r 

Iri— rl r 



(B2) 



(B3) 



(B4) 



-it'{Ho-E) 



ri=0) . 
(B5) 
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By introducing the Heisenberg picture operator 



ri(t) = e*^o*rie-^^°* , 
Lp(r) may be simply expressed as 

r 



(B6) 



1 



X 



|ri(t') 



|ri=0) 



(B7) 



where we have used (p| {Hq ~ E) = 0. Since the Hamihonian Hq describes a free particle, 
the equation of motion for ri{t) gives 



Pi. 



ri(t) =ri(0) + ^^t. 



m 



(B8) 



We shall need only the leading order term of Lp{r) for p being much larger that l/ao- 
This is the term that gives the leading correction of order rj'^ = l/(]9ao)^. It is then valid 
to treat the momentum operator pi as commuting with the coordinate operator ri since 
[ri, pi] ~ 1 ^ paQ. Thus, the leading order in rj evaluation of the matrix element in Eq. (B7) 
is given by its classical limit with the operator pi replaced by its eigenvalue p and ri set to 
zero and with (p|ri=0) = l: 



^p(r) 



oo roo 



1 



dt dt' { — , 
^0 \\l{tM')-r\ 



1 



At; 



-{t+t') 
m 



1 



I m I 



(B9) 



Performing a change of variables 



mr , „ , t' 

s=—{t + t'), A = — -, 
p t + t' 



(BIO) 



puts this in the form 

,2 



^p(r) 



•fyi'^ roo ( \ 1 

/ dss / dX\— — rA-ufsr) 

p^ Jo Jo \\ps — r| s 



1 



1 . 



(Bll) 



|pAs — r| 

Since only the angular average in r contributes, it is convenient to use the spherical 
harmonic expansion 
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oo I 



1 r: 



(B12) 



to obtain 



^p(r) 



oo I 



X 



oo I' 



An s' ' 

X! rfj/ I 1 ,if+i^*m'(p)^i 



Af (sr) 
(B13) 



where 



s< = mm{s, 1}, s> = max{s, 1} , 
s'^ = min{As, 1}, s'-> = max{As, 1} . 

The orthonormahty of the spherical harmonics and the property 



(B14) 



21 + 



(B15) 



now yield the angular average 



-— ^ ^2 ^oo H 

Lp(r) = / dss / dA 

Jo Jo 



(B16) 



For the first term in the integrand, exchanging the order of the integrals and sum, it is 
straight forward to evaluate 

J r.ll 



/■I °° 1 1"^ S S 

Jo ^^g2lTTio 

Jo ^^21 + 1 \Jo 

„i oo / \/— 1 

= / ^AE^' 

•^0 1=1 



21 21 + 2 



2 2 



(B17) 



Therefore, 



^p(r) 



m 



X /"OO 

2 



/oo 
(is (In s — s + 1) r A'u(sr) 



(B18) 



where we have carried out the remaining A integral. The expansion (|B12| ) gives 
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Av{r) = {dh')pf,{r') 



1 1 



(B19) 



Inserting this expression for Av{r) into Eq. ( [B18| ) and performing the s integral produces 
,2 r 



^p(r) 



m 



2p2 



{dh')pf,{r')e{r'-r) 



(B20) 



Putting this result in to Eq. ( [BBD gives 

Jf^ = -ll- f{dh)pf,{r) hdh')pf,{r') 



^-llnM^ 
2r' 4 \r 



(B21) 



where the fact that the integrand is symmetric in r and r' has been used to replace 9{r—r') 
by 1/2. 

To simplify this result, we note that 



r 



to obtain 



ao 



mU\rH,-Hfr\i) 
:^l/)(/|kp^l |^) + (^|/)(/|^) 



Since (r|z) ~ exp{— r/ao}, 



ao(2|/)(/|(2r-p + -)K) + (2|/)(/K). 



ia^i ■ p|i) = — |i) 



and so 



(^l^l/)(/|rK) = (^|/)(/|^K) 



(B22) 



(B23) 



(B24) 



(B25) 



or 



j {d^ri)pfi{ri) J {dh2)pfi{r2) (^^^ = almVfi 



(B26) 
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Therefore, Jfi can be repackaged as 



with 



(B27) 



(B28) 



This is the result (|3.2(]| ) quoted in the text. 

We shall calculate Kfi for the final ^He"*" state being in Is, 2s, and 3s states. To evaluate 
Kfi, it is convenient to consider the integral 



/(e) = / {dh,)pf,{n) J {dh,)pf,{r,) {^fj , 



(B29) 



which enables us to express Kfi in terms of /(e) as 



1 



4 de^ 



(B30) 



e=0 



For the case where the final ^He^ is in its ground state, the Is state, using the wave functions 



1 / 2 \ 



1 _ 



an e 



7"o 



(B31) 



we have Pfi{r) oc (f)*j{r)(j)i{r) oc e with s = 3/ao. Using this form of Pfi{r) and performing 
the changes of variables Xi = sri and X2 = sr2 yields 

1 poo POO 1 

/(e) = / dxixl-'e"^ / rfx2X^+^e-^2 = -1(3 + e)r(3 - e) , (B32) 
4 JO Jo 4 

where we have used /(O) = 1 to get the correct normalization factor. Using 



r(l + e)r(l-e) 



7re 



smvre 



(B33) 



expression ( [B30[ ) now reads 



~ 1 d^ r7re(4-e2)(l-e2) 

Kli = — -tt: 



16 de^ 



sm vre 



e=0 



V _ 5^ 
12 ~ 8 



(B34) 



Similarly, exploiting the wave functions for the final 2s and 3s ion states 
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02s(r) 
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1 \^/^ 



1 /2_\3/2 



r 

1 



81v37r Vao/ 
we find the corresponding and K^i. 

I 12 8/ ' 



27 - 36 ( — + 8 — 

cto/ Vao 



exp 
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Bar 



TT 



i^3. = I -T 



1045\ 



12 648 ; 



(B35) 
(B36) 
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2. Evaluation of T? for final atomic states with non-zero angular momentum 



Here we provide the intermediate steps for deriving the result ( ^.4|) from the expres- 
sion (|5.2| ) for T^. Following steps similar to those that lead to Eq. ( P 1 1|) in the previous 
subsection, we can rewrite to the leading order as 



//■oo 1 
{dh2)<P}{r2)Ur2) ds- — . 
JO |ps — 

For the final state (/| with angular momentum (/, m) and thus the wave function 



r) = R},{r)Yav) , 



we use the expansion (|B12|) to get 



I'm' 



-I'm dr2rlR}iir2)Mr2)7^r-TyLiP) ds 6(1-3)3' 
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21 + 1 

roo 



(3-1)3 



-i-i 
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APPENDIX C: INCLUSIVE DETAILS 



1. Coulomb Wave Functions 



To facilitate the calculations below, we display here the Coulomb wave function for the 
scattering state on the final ^He''"'' nucleus of charge two |^ 
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iJk{r) = (r|k) = e"^r(l - 2ij)e^^-''F{2i-f, 1; ikr -ik-r), (CI) 
where F{a, c; z) is the confluent hypergeometric function, and the parameter 7 is defined by 

^ kaQ ^ ^ 

The appropriate integral representation of the confluent hypergeometric function yields 

Jc 2Tii 

where the contour C wraps the cut in t plane connecting the branch points t = and t = 1 

in a counter-clockwise sense. This wave function has the partial wave expansion 

00 

Mr) = E(2^ + l)Pi{cos9)e'^'/^Rki{r) , (C4) 

(=0 

with 

Rki{r) = ^^'^^^^^~^^'^\ 2krye''^^F{l + 1 - 2^7, 21 + 2; -2^kr) . (C5) 

Here Pi(cos^^) is the Legendre function of the first kind and 6 is the angle between k and r. 
These radial wave functions can be recast into the following integral forms by using different 
integral representations of the confluent hypergeometric functions 



ft,(r)^e-r(l + /-2n)(35^£|-e-'(*+5) {t - -) , (C6) 

and 

Rkiir) = -e-"T{-l - 2z^){2kry £ —e'^'^' (t + - j (t - -j , (C7) 

which we shall use later. Here the contour C wraps the cut connecting the branch points at 
t = ±1/2 in the counter-clockwise sense. 



2. (k|z) 

The matrix element (k|z) can be evaluated by using the Coulomb wave function directly. 
Since the wave function of \i) is spherically symmetric, in view of Eq. ( |C4[ ) only Rko{r) is 
needed. Utilizing Eq. ( |C6|) for / = and performing the r integral yields 
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(^Ik) 



Anr'^dr (f)i{r) Rko{r) 



_3/2e"^r(l-2i7) / dt 1 



t + 



-1+2*7 / -l-2i7 
t - 



(C8) 



{2ik)^ Jc 27Ti (t+i7/2)2 

Deforming the contour C to wrap the double pole at t = —i'y/2 enables us to evaluate the 
contour integral as 



|k) = —A^/^la, 



3/2 e-^r(l-227) d 



Consequently, 



(m 



3/2 



e"^r(l-2i7)- 



1 + 7')' U7 + 1 



2i7 



t=-i'y/2 



(C9) 



(2|k) 



2567r2 1 76 
l_e-4^7fc3 (1 + ^2)4' 



-87 cot ^ 7 



(CIO) 



3. K"{0) 



In this subsection, we compute 



K"{0) = f-{t\r-%){k\r'\i) 
de^ 



(Cll) 



e=0 



for Coulomb states in the continuum |k). We first express fC"(0) in terms of several one 
parameter integrals, then study its asymptotic behavior as a function of 7, and finally present 
numerical results. 

Using the representation 

and the Coulomb wave function (|C6D for / = and performing the radial integral, we get 



(«i I ik> 



7ra§e''r(l - 2i',) .~ 

z : 7 / dXX 

2r 1 + e ' Jo 



X 



dt 
c 2m 



t + 



27(1 + x) 



l \ -l+2i7 / -1-2*7 

'+2 '-2, 



(C13) 



Deforming the contour to wrap the triple pole at t = —^7(1 + x)/2 gives 
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2r(l + e 

1 







which can be also written in form 



' + 2) ('-2 



(C14) 



t=-i7(l+x)/2 



(^1 



|k) 



r(l+e) ^ io 



1 



[l+72(l+x) 



-47 cot ^ 7(1+0;) 



. (C15) 



Putting this into the definition of K(t) yields 

o sin vre 47r7 4 /""^ , /■°° , / x 
vre 1 - e-^'^T 7o 7o V?/ 

rf2 



K[e) = Wnal 



[1+72(1+0;) 



-47 cot ^'y{l+x) 



1 



^—47 cot ^-f{l+y) 



(C16) 



Rescaling the dummy variables x,y hj a. factor 7 and taking the derivative with respect to 



e gives 



K"{0) = -32TTaf,- 



2 ^,2 



7 



(C17) 



where we have defined the integrals 



POO 



-47 cot -'-(7+x) 



(C18) 



_1 + (7 + 

with k = 1,2. We have exhausted our analytic power at this point. This expression is the 
starting point of the numerical computation. 

Before displaying the numerical data of K"{0), let us explore its asymptotic behavior as 
7 — > and 7 00. For a small 7, we first study the corresponding asymptotic behavior of 
Ik- Noting that 



cot ^(7 + x) = — — tan "^(7 + x) 



(C19) 



expanding the integrand, and keeping only the first two leading terms in the small 7 expan- 
sion, reduces Eq. (|C18|) to 
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4(7) 



ax^ 



1 



27a; 



, 47 

+ tan X 



(C20) 



Examining Eq. ( |(J17| ), we find that, to the first two leading terms in the small 7 asymptotic 
expansion, we need to keep only the first term for 12(7) in the equation above. For /2, the 
leading term can be calculated by doing the x integral by parts and evaluating a contour 
integral. The result is 



-TT 



as 7 ^ - 



(C21) 



For /i, the two leading terms must be kept. It can be expressed as 

1 



/i(7) ~ 2e-2"^ / dx 
Jo 



+ 



47 , _i ^, 27X 
tan x + mx- 



(C22) 



;i + x2)2 (l + x2)2 dx^l+xY_ 
where we have split the integrand into two parts and done the x integral by parts for one 
of them. The x integrals in the equation above can be evaluated either by using their 
corresponding indefinite integral or by writing them as derivatives of the Gamma functions. 
We shall not go through the details of the derivation here, but only display the results: 



TT 



TT 



+ 37 as 7 — * 



(C23) 



2 V 2 

Inserting these into Eq. ( |C17|) and making appropriate expansions of the other functions of 
7, we obtain the asymptotic form of K"{0): 



K"{0) ^ -87r=^7^a|^ + 1607r^7^a;; as 7 ^ . 



2, ,7^3 



(C24) 



We now study the asymptotic expansion of K"{0) for large 7. It is convenient to look at 
the 7 00 limit of Eq. ( |Cl(j| ), which is 



K(e) = 647rVo^^ H dx H dy ( - 
ne Jo Jo \v, 



X ■ 



d^ 

dx"^ 



1+x 



-4/(l+x) 



rf2 



-4/(1+!/) 



dy^ \{i + yf 



(C25) 



Taking the derivative with respect to e produces 

"27r2 



K"{{]) = -1287rVo 



-8 



(C26) 
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where Ci,C2 are two constants 

Numerically, Ci = 0.0525 and C2 = —0.0664. Using these values, the asymptotic behavior 
becomes 

K"{0) -3.19 -fal as7^cx). (C28) 




FIG. 11.: Numerical curve for K"{0)/'y as a function of 7. The dashed horizontal line is 
its asymptotic behavior for large 7. The unit for the vertical axis is a^^ and the horizontal 
variable 7 is dimensionless. 

We present the curves of K"{0) as a function of 7 and its large 7 asymptotic behavior 
in Fig. 11. Since the small 7 asymptotic form ( P28|) fits K"{0) only at a small region, we 
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FIG. 12.: Numerical curves for K"[0)/^ as a function of 7 and its asymptotic behavior. The 
lower curve represents K"{0) and the upper dashed curve represents its asymptotic behavior 
for small 7. The unit for the vertical axis is Qq^ and the horizontal axis represents the 
dimensionless variable 7. 
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FIG. 13.: Relative accuracy of the asymptotic form [K"{0) + 877^7^ - 1607r27^]/(87r376) for 
small 7 as a function of 7; both vertical and horizontal axes are dimensionless. 
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display it separately in Fig. 12 and Fig. 13. Our numerical data shows that at the large 
7 region, K"{0) quickly approaches its asymptotic form ( p28| ). For 7 larger than 10, the 
asymptotic form has the accuracy of 0.01. At the small 7 region, one can see from Fig. 13 
that the asymptotic form ( p24| ) has the accuracy 0.01 for 7 smaller than 0.02. 



4. i^2b 



We now turn to evaluate 



K2^ = Y.iE^-EfrK'\0) 

b 



(C29) 



e=0 



where the sum runs only over the final bound states. This was defined by Eq. ( |4.26|) in the 
text. We first note that 



{E,-Ef){f\r^\z) = {f\{r^H,-Hfr^) \z) 



{f\r' 



/' 

e' , 6(e + l) 



H ir ■ p 



mr 



(C30) 



Since 0i(r) oc exp(— r/ao), the operator ir ■ p is equivalent to the factor — 1/ao- Therefore, 



iE,-Ef){fy\t) = Ify{f\r^ 



2(l-e)ao e(6 + l)ag 



(C31) 



Inserting Eq. ( P31|) into the definition (|C29|) and carrying out the derivatives is straightfor- 
ward algebra. The result may be simplified by two remarks. Since the radial Schrodinger 
equations are real, so are their regular solutions, except for an overall phase factor. Hence, 



mr)\f){f\9{rm = mr)\f){f\f{rM . 



(C32) 



Moreover, the definition ( p29D is not changed by the reflection e ^ — e . Making use of these 
remarks puts the result in the form 



Y.{E.-EjrK"{0) = VE |-2(^|4|/)(/|4n) + m-\f){f4\^) - 8{^\^\f){f\^\^) 
b fel^'^ r r J. J. 
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d , , / r 
+8- (z| - 
ae \ao 



-l-e 



-2+e 



e=0 



-l+e 



(C33) 



e=0- 



The first three terms in the right hand side of Eq. ( |(J33| ) may be readily calculated by 
using the bound state wave functions for the Coulomb potential. The bound states for 
Coulomb potential are labeled by the principal quantum number n and the angular quantum 
numbers l,m. Since only the final S-wave states give non-vanishing matrix elements, we can 
simply use the principal quantum number n = 1, 2, ■ ■ ■ to denote the final states. The wave 
function may be expressed in terms of the generalized Laguerre polynomials L'^{x) which 
may be written in terms of an integral representation to express 

1 r dt / 1\" 

^i2^"''p^"^^+^^^^^'"/''''°^H^+i) ' ^^^^^ 

with the contour C encircling t = in the counter-clockwise sense. Using this and the initial 



1 / 2 \ 
^ ' ' y/n \aoJ 



wave function ( P31|) and performing the trivial radial integral gives 



ao 



I/) 



2^/2 r(3-m) 



c 2T(i ( ~'~ t ^ 



2 4t\"-3 
- + - 
n n J 



(C35) 



The contour can be deformed to wrap the pole at t = — n/4 — 1/2 to carry out the t integral. 
For m = 2, there is a contribution from the contour at the infinity. Therefore, 



ao 
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2^/2 1(3- 



-m) /n 
^5/2 (2_^ \1 



3—m 
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For m = 1, 2, it is straight forward to evaluate 



(^1^1/) 
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8V2n 



n 
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(^1(7) I/) 
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2s 3/2 

nj 
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n + 2 
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m,2 



(C36) 



1 



n - 2 
n + 2 



(n + 2)2 \n + 2 

Hence, the first three sums in the right hand side of Eq. (|C33| ) can be written as 

nn 2 



(C37) 
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E(^l§l/)(/l%) = 8E:^ 



n=l 

00 



1 - 



On,., „^ 1 
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n + 2 
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{ n{n + 2)2 \n + 2 

n /n-2\2("-i) 
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V ^ ^ ^i{n + 2Y\n + 2 



n-2 

n + 2 

1.58. 



4.78 



(C38) 
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For the fourth and fifth terms in Eq. ( |C33|) , we shall compute the bound-state sums by 
exploiting the completeness of the the sum over all the intermediate S-wave states. We write 
the bound-state sum in terms of the identity operator (in the S-wave sector) — which is 
equivalent to using the closure result — and then subtract the sum over all the continuum 
S'-wave Coulomb scattering states: 



El/)(/l = i 



1 



27r2 



k^dk\k){k\ . 



(C39) 



Since the closure results — the terms involving the identity operator — do not have any e 
dependence and thus do not contribute to the derivative, the identity operator in Eq. ( |C39|) 
may be omitted in our calculations. 
Thus the fourth term in Eq. 

l+e 



2-€ 



e=0 



may be written as 

(.|(^)"f.)(.|(^)^-fo 



°^dk—- 
27r2rf6 



e=0 



Using the representation ( |C12| ) and taking the derivative with respect to e gives 



(^lf^)"l/)(/lf^)"'l^') 
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4717 



4 f°° 2 
n Jo l-e-4^T 



J7 g-47Cot-7J3(^) + (l-e-4^™*-^^)/4(7) 

1+7'' 



where the integrals I3 and J4 are defined by 



and 



h{l)= dx\nx— — ^— — 
Jo l + 7^(l + x)"' 



hiri) = / dx\nx 

Jo \l + Y[l 



-47 cot 7(l+x) 



^—47 cot 7(l+x) 



[1+72(1 + X)2]2 

We display here the numerical value of this fourth term 



(C40) 



(C41) 



(C42) 



(C43) 
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In the same fashion, we can deal with the fifth term in Eq. (|C3 



e=0 
oo ][ 
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Evaluating the derivative gives 



TT^ /-oo 



X 



3 Jo 27r2 
7^(l-x) 



(.1^1^) 
r 



e=0 

/•CO /-oo / 7- \ i-CO 1 ZtTTO/ 



-647ran / (ix / c/yln — , , 



_ — 47Cot 7(1+0;)— 47 cot 7(1+?/) 
[l+72(l+x)2]2[l+72(l+^)2]2 



8 r°° , 47r7 
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where we have defined the integral 
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and used the result 
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/O [1 + 72(1+X)2]2 

Numerically, this fifth term has the value 



2(1+72) 



-47 cot 7 
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2.02. 
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Putting these numerical values into Eq. (|C33| ) yields finally 



(C49) 



K2^, = Y.{E, - EffK"{0) ^ 32.26%' . 
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5. M(k) 



We turn now to calculate M{k) defined in Eq. ( p.8|) , 

M{k) = {E, - E,f f drT'RUr)<i^,{T] 



(C51) 



We shall express M{k) as a one parameter integral. 

To do this, let us consider first the simpler sum S,") defined by 

oo 

G(e',n = E(2/+ik.Kmz(n, (C52) 
1=1 

where 

roo 

g,,[^)^k-^ drRUr)<P,{^,r) , (C53) 
Jo 



4>,(^,r) = ^e-^^/^\ (C54) 



with 



which, up to an overall constant, is the ground state wave function of a hydrogen-like atom 
with nuclear charge +^e. The reason for considering G{C,',C,") is that, as we shall see later, 
M{k) may be expressed in terms of G{^',^") and G{C,',C') i^^iy be evaluated in a closed 
form. We first complete the evaluation of G{^',^"). Using the orthogonality relation of the 
Legendre function. 



/' dxPi{x)P,{x) = 777^5/,/' , (C55) 
J-i 2< + i 

the expansion ( P^ ) expresses the radial Coulomb wave function as 

Rki{r) = A /' rfcos^Pz(cos^)V^k(r,^,0). (C56) 

^ 1 J — 1 

Inserting this into the definition of QkiiO ^"^^ using the completeness relation 

oo 

$3(2/ + l)Pi{x)Pi{y) = 2S{x - y) (C57) 

1=0 

yields 
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where the second term accounts for the fact that we do not include the / = term in 
the summation ( |C52| ). Making use of the integral representation ( |C3| ) and the definition 
of 0i(^,T), it is straight forward to first perform the r integration and then, deforming the 
contour of the t-integration to encircle the pole as in the evaluation of (|C9D , to compute 
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The remaining angular integration in Eq. ( P58|) may now be done in closed form: 
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. For the second term, 



2(e' + e 

With this, it is straight forward to evaluate the first term in Eq. 
using the integral representation ( |C7| ) for the radial wave function Ruir) for Z = and doing 
the r integral, we find that 

g^7r( 1-2^7) r dt 1 
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where the contour integral has been evaluated in the fashion as in the previous evaluations 
in Eqs. (|C9|). Note that, while deforming the contour, the contour at infinity contributes the 
1 in the last parentheses. We have now found that 
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As mentioned before, G{C,',C,") is related with M{k), as we shall now prove. To facilitate 
the derivation, it is convenient to define 
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which generates the desired matrix element involved in the definition (|C51|) via 
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The Schrodinger equation satisfied by Rki{r 
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enables us to relate guiCj and fkiiCj through 
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where we have done the integral by parts and used the explicit form of the wave function 
(f)i{C,,r). Relation ( p66| ) may be easily inverted: 
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where the boundary condition has been appropriately chosen. Combining this with Eq. ( |C64|) 
yields 
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Consequently, M{k) may be expressed in terms of G{^',^") as 
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(1+7^)2 

where we have used the fact that G{^',C,") is symmetric about its two arguments. 

Inserting the exphcit form (|C62|) for G{^',^") into the expression above for M{k) and 
making the changes of the variables ^ = 1/t, ^' = 1/t', and C," = 1/t" yields the lengthy 
expression 



■47 tan 7 



1 + 72 '^'^ ' (1 + 72)2 

where the integrals 6*1(7) and ^2(7) are defined by 
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The integral defining (^2(7) may be simplified to a one parameter integral. This can be done 
by first making the change of variable t' = t, t" = st, which gives 
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where we have done the s integral by parts for the second term. We then do the t integral 
by parts and observe that the s integral may be completed in a closed form to obtain 
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t Jo s \[l + (t/7)2][l + (st/7)2] 
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1 
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[1 + (tsh) 
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+- 



1 dt 



^47tan-i(V7) 



A Jo 



-l)(e 



47tan-l(i/7) 



-8 ^ ?J 

1 / l_\ dt 



47 tan 1 7 ^ _ ]^ N _|_ 2 



- 1 - 4t) 

1 dt 1 

Ti + t 



^g47tan l(t/7) _ -j^^ 



^47tan-i(i/7) /'p47tan-l 7"^ 
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where we have chosen t as the one parameter integrals variable uniformly and have appro- 
priately performed the t integral by parts several times. Thus, M{k) can be written in terms 
of one parameter integrals: 



I67I 



^- -47tan"i7"i 
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lo t{t+l) 

We now study the asymptotic behavior of M{k). For small 7, we shall only keep the 
terms up to order 0(7'^). This may be done by using the following asymptotic expansions of 
the two integrals involved in Eq. (|C75D : 



^ _^^_g47tan-i(t/7)-47tan-i7-l ^ ^ 



(l + t)2 

1 dt 

t{l + t) 



^47tan-i(t/7) 



- + 47^ In 27 
: — 27r7 In 27 . 



(C76) 

(C77) 
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These results can be derived by doing the t integral by parts and expanding the exponential 
into a power series. Putting them into the expression for M{k) gives 



M(fc) ~ 27r7-(4+7r^)7^+167r7^1n7+(20+161n2)7r7^ % X as 7 -> . (C78) 

For 7 being large, it is not hard to see from the expression (|C75|) that the leading behavior 
of M{k) is linear in 7. The proportion constant is evaluated numerically to be 0.360. Thus, 

M{k) ~ O.36O7 as 7 ^ 00 . (C79) 

We now display the numerical curves of M{k)/'y and its asymptotic behavior as functions 
of 7 in Fig. 14 and Fig. 15. 

APPENDIX D: DETAILS OF THE EXCHANGE TERMS 

We first examine the third term T3 in Eq. ( |2.27] ). For this term T3, only the "vacuum 
state" |0) (i.e. the tritium nucleus) contributes to the intermediate states when inserting a 
complete set of states just before the operator A''p(0). Hence, 

= -^}(0)^(p|{^ + .(r))|.>. (Dl) 

where (p| and \i) are the electron states of the final beta ray and the initial tritium atom. 
The Schrodinger equations for these energy eigenstates can be used to rewrite 

T3 = -0}(O)(pK). (D2) 

This is an exchange term which represents the amplitude for the electron produced by the 
weak interaction shaking the bound electron of the original tritium out and being bound by 
the final helium nucleus itself. We shall not evaluate T3 because it is canceled by a piece of 
the exchange term which we now turn to discuss.0 



^^It is not hard to see that T3 is of order t]^ by examining the expression ( |D1| ) written in terms of 
wave functions. Since, for any reasonable comparison potential, the integral 
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FIG. 14.: Numerical curve for M{k)/'~f as a function of 7. The long-dashed horizontal line is 
the asymptotic limit for large 7 given in Eq. ( |U79| ). The short-dashed curve at small value 
of 7 is the approximation given in Eq. ( |(J78| ). The vertical axis is in units of Ry^c^o while the 
horizontal axis is dimensionless. 
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FIG. 15.: Numerical curve for M{k)/'y (solid curve) as a function of 7 compared to the small 
7 expansion (dashed curve) given in Eq. ( |C78| ). The vertical axis has the unit ^/^al and 
horizontal axis is dimensionless. 
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The exchange part in T2 of Eq. ( p.5|) involves the ket \i, r2=0) For the evaluation of 
this exchange term, it is convenient to interchange the roles of ri and r2 in Eqs. ( p.7|) , (|3.8| ) 
and (13.91) and define 



Hi 
Hi 



P? e2 



2m Anri 
2m 47rr2 

„2 



(D3) 
(D4) 
(D5) 



47r|ri — Attti 

This enables us to expand Tg in powers of Hj as was done before [c.f. Eq. ( |3.10|) ]. The 
leading term in Tg, which we now denote as T^, is 



2e' 



47r|ri — r2| 47rri 



i;(ri) 



H2 + Ei- E - Ef -le 



-|r2=0) 



(D6) 



upon using the fact that \i) is an eigenstate of Hi. Noting that (/| is an eigenstate of H2, 
we may write 



(rf3ri)C(ri)0.(ri)(/| 



^47r|ri-r2| 47rri / H2+Ei-E-Ef- 



-le 



|r2=0) 



+ v{ri) 



(D7) 



The second term here exactly cancels T3 in view of Eq. (|D1|) . 



converges, and 1/p is much less than the Bohr radius, we can replace the wave function (f>i{r) by its 
value at r = to obtain the leading order contribution. Noting that 0j(O) and 0i(O) are of order 
(am)^/^, and approximating the integral displayed above by the Fourier transform of the Coulomb 
potential evaluated at momentum p, yields 



T3 ~ (am)3/2 (eVp^) {amf/'^ ~ rj^ . 
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To evaluate the sum T^^ + T3, the first hue in Eq. ( P7|) , we follow the method of Appendix 
B in writing the denominator as the integral of an exponential, use the Heisenberg equation 
of motion with {f\{H2 — Ej) = 0, and express the Coulomb Green's functions as Fourier 
integrals, to obtain 

In leading order, which is equivalent to the p —>■ 00 limit, we may replace the outgoing beta 
electron wave function 0p(ri) by the free plane wave e~^^''^^, neglect Ei relative to E = p^/2m, 
and replace r2{t) by the free particle motion, 

r2(t) = rs + pt/m. (D9) 

Ordering the resulting exponential gives 

g-ik-r2(i) ^ g-ifc2t/2mg-itk p/mg-jk r2 (DIO) 

Since the momentum operator generates a spatial translation, we now arrive at 

poo 

J rfte'P'*/'™ {e"*^'*/''"(/|r2 = kt/m) - (/|r2=0)} . (Dll) 
To exhibit the leading order contribution, we change variables by writing 

Tn 

k=-r, (D12) 



and 



TTIT 1 

t= r, ri = r + -u. (D13) 

p p 



This change produces 

J Jo T 



^iprr/2 |e-w/2r^*(r) - 0}(O)} j + u/p)e-'^P-^/^>^ . (DM) 



In the p — s> oo limit, 0i(r + u/p) may be replaced by 0i(r), and we encounter 



/|^e-(P-^/-)- = 5(^)(p-r/r) 



= 6il-T)J2Y;^ip)YU^), (D15) 

Im 

which produces 

T° + T3 ^Wal f {dh)\ Y: Y;^,{m'm'{^) {0}(r) - e'^'-0}(O)} 0.(r) . (D16) 

•' ^ I'm' 

The last term involving e*^^0j(O) vanishes in the large p limit by virtue of its infinitely rapid 
phase oscillation. The angular part of the r integration just picks out the / and m values of 
the final atomic wave function, 

0}(r) = i?},(r)l^:„(r). (D17) 
Accordingly, to the leading order, the exchange amplitude is given by 

T'^+T,c^iATxr^'alYai>) drR},{r)Ur) , (D18) 

which is of order rj^. 



Although the amplitude (|D18|) is of order 77^, it contributes to the decay rate only through 



the order 77^. For 1^0, this result is of order rj'^ relative to the leading direct result (|5.4| ) which 
is of order rj. Therefore the exchange effect gives a correction to the decay rate at order t]'^. 
For case 1 = 0, the leading exchange amplitude (P18|) is relatively imaginary. Since the 
leading imaginary amplitude appears at the order rj, the exchange correction to the decay 
rate is only of order 7]^ We note here that this contradicts with the results in the literature 
14|, where the dominant exchange amplitudes appear in order 77^. 
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